On the effect of the thermal gas component to the stabihty of vortices in trapped 

Bose-Einstein condensates 
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We study the stability of vortices in trapped single- 
component Bose-Einstein condensates within self-consistent 
mean-field theories — especially we consider the Hartree-Fock- 
Bogoliubov-Popov theory and its recently proposed gapless 
extensions. It is shown that for sufficiently repulsively inter- 
acting systems the anomalous negative-energy modes related 
to vortex instabilities are lifted to positive energies due to par- 
tial filling of the vortex core with noncondensed gas. Such a 
behavior implies that within these theories the vortex states 
are eventually stable against transfer of condensate matter 
to the anomalous core modes. This self-stabilization of vor- 
tices, shown to occur under very general circumstances, is 
contrasted to the predictions of the non-self-consistent Bo- 
goliubov approximation, which is known to exhibit anoma- 
lous modes for all vortex configurations and thus implying 
instability of these states. In addition, the shortcomings of 
these approximations in describing the properties of vortices 
are analysed, and the need of a self-consistent theory taking 
properly into account the coupled dynamics of the condensate 
and the noncondensate atoms is emphasized. 



PACS number(s): 03.75.Fi, 05.30. Jp, 67.40. Db, 67.40.Vs 



I. INTRODUCTION 

The realization of Bose-Einstein condensation in 
trapped, dilute atomic gases has opened up unique 
ways to investigate how weak particle interactions affect 
the quantum-statistical phenomenon of bosonic conden- 
sation. Especially interesting is the question whether the 
interactions can sustain superfluidity in such systems. A 
characteristic property of superfluids is the ability to sup- 
port dissipationless flow, a feature which is manifested 
in microscopic length scales as the stability of quantized 
vortices even in fluids confined by nonrotating vessels. 
This aspect on the superfluid properties of dilute atomic 
Bose-Einstein condensates (BECs) has been under vig- 
orous theoretical investigation during recent years ]^|2^] 
(for a review, see Ref. p9|), and has been further moti- 
vated by the experimental advances in creating and ob- 
serving vortex structures and their dynamics in these sys- 
tems |po|--p6|. However, the majority of the analysis has 
been carried out within the zero-temperature Bogoliubov 
approximation, which neglects the effects of the noncon- 
densed gas component always present in interacting sys- 
tems. Especially, in addition to providing a dissipative 



mechanism for the condensate vortex state ^] , the non- 
condensate component also modifies the excitation spec- 
trum of the system, and may thus in principle affect the 
stability of vortices. 

In discussing the energetic stability, one has to distin- 
guish between global, thermodynamic stability and sta- 
bility against small deformations, or local stability. A 
physical system is defined to be locally energetically sta- 
ble if its state is a local minimum of the free-energy func- 
tional under the given external constraints — if the state 
is also a global minimum of energy, it is thermodynami- 
cally stable [Q. In studying the superfluid properties of 
dilute atomic BECs, it is crucial to find out whether these 
systems support persistent currents without application 
of any external drives — especially, whether vortices are 
stable in condensates confined by nonrotating, station- 
ary traps. In this context, the relevant stability criterion 
is the local energetic stability, because such vortex states 
are not thermodynamically stable even in superfluids — 
the global energy minimum in a stationary trap always 
corresponds to a nonrotating state. On the other hand, 
local stability of bosonic systems is associated with the 
positivity of the energies of the normal modes. Our inter- 
est hence is to determine the local stability, or metasta- 
bility, of vortex states by investigating their normal mode 
excitation spectra. 

The local stability of vortices in harmonically trapped 
atomic BECs has been addressed mainly using the zero- 
temperature Bogoliubov scheme, in which the noncon- 
densed gas is neglected. It has been shown that within 
the Bogoliubov approximation (BA), vortex states in 
nonrotating traps always have anomalous excitations 
with positive norm but negative energy (with respect to 
the condensate state), which essentially correspond to a 
translational displacement of the vortex from the cen- 
ter of the condensate and/or bending of the vortex line 
|6[-^, p^|j2C|j2^ . The negative energies of these excitations 
imply that vortex states are energetically unstable: In the 
presence of any initial perturbation and dissipative mech- 
anism, these modes would start to be excited through 
the transfer of matter from the condensate state. In ad- 
dition to possible bending and deformation of the vor- 
tex, the negative-energy displacement instability mode 
implies that the vortex is prone to start spiralling away 
from the center of the condensate Q . By computing the 
energy of the condensate particles as a function of the 
vortex displacement from the condensate center it 
has been verified that it is energetically favourable for the 
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vortex to eventually spiral out of the condensate and be 
annihilated. In conclusion, when the thermal gas com- 
ponent in the system is neglected, the mean-field theory 
predicts vortices to be energetically unstable. 

At temperatures comparable to the condensation tem- 
perature Tbec, a substantial portion of the gas is non- 
condensed and the validity of the instability prediction 
given by the BA is not obvious. However, surprisingly 
enough, for vortex states its validity is questionable also 
in the zero-temperature limit. In interacting systems the 
noncondensate density never vanishes exactly, but the 
residual noncondensate density in the zero-temperature 
limit is practically negligibly small for weakly interacting 
condensates. However, because of the existence of the 
negative-energy modes within the BA for vortex states, 
in the presence of dissipation a nonnegligible portion of 
the condensate always "leaks" to noncondensate popu- 
lations of these anomalous modes. The consequent non- 
condensate density concentrated in the vortex core re- 
gion, in turn, may in principle affect the local stability 
of the vortex state by changing the anomalous mode en- 
ergies themselves. A similar phenomenon is observed for 
two-component condensate systems, in which an irrota- 
tional condensate component can stabilize vortex states 
p^ . The crucial property determining whether such a 
stabilization is also possible due to finite noncondensate 
density, is the amount of noncondensate required to lift 
the anomalous modes to positive energies: If the lifting 
effect is too weak, the whole condensate vortex state will 
leak to macroscopic occupations of the anomalous modes, 
and the vortex will be destroyed. 

Previously, the self-stabilization mechanism described 
above has been demonstrated in the case of an infinitely 
long, axisymmetric, static vortex line located in the cen- 
ter of a harmonically trapped condensate. It has been 
shown that within the self-consistent Popov approxima- 
tion (PA) and its recently proposed gapless exten- 
sions Gl and G2 |^,^, the noncondensate in the vor- 
tex core lifts the anomalous modes to positive energies 
for thermally equilibrated vortex states even in the zero- 
temperature limit ||l^,^. However, even slightly off- 
axis vortices are not stationary, but rather move about 
the condensate typically at velocities, for which the as- 
sumption of adiabaticity implicitly made when applying 
stationary self-consistent mean-field theories to describe 
time-dependent phenomena is highly questionable [ p7| . 
This raises the question of applicability of the results 
obtained for axisymmetric vortex states to the general 
off-axis vortex case. In this paper we generalize the pre- 
vious results for the PA, Gl and G2 by showing that 
the anomalous modes of a circularly precessing off-axis 
vortex line penetrating a sufficiently strongly interact- 
ing, trapped condensate are lifted to positive energies as 
the noncondensate partially fills the vortex core region — 
the noncondensate fraction required for such stabilization 
depends on the effective strength of the repulsive inter- 
actions in the system, but it is of the order of the ratio 
of the vortex core volume to the volume of the whole 



condensate. 

In conclusion, the Bogoliubov approximation predicts 
vortices in harmonically trapped BECs to be energeti- 
cally unstable, while the self-consistent Popov, Gl and 
G2 mean-field approximations (and other approxima- 
tions having a similar structure; we refer to them with the 
term Popov-like approximations (PL As)) predict them to 
be stable. This discrepancy can be contrasted to the ex- 
perimental observations of vortex precession in spherical 
condensates Theoretically, it has been shown that 

the precession direction of a vortex is related to the en- 
ergy of the vortex displacement mode: If the mode is 
anomalous with a negative energy, the vortex precesses 
in the direction of the condensate flow, and vice versa 
po[-^. Except for a minority of so-called rogue vor- 
tices, the precession direction observed supports the pre- 
dictions of the BA rather than those of the self-consistent 
approximations. As the results derived in this paper sug- 
gest that the non-adiabaticity in the time-development of 
moving vortices is not responsible for the discrepancy be- 
tween the experimental observations and the predictions 
of the PLAs, we are essentially left with two possibil- 
ities: Either dissipation in the experimental conditions 
has been too weak for the vortex states to be sufhciently 
properly thermalized, or these self-consistent approxima- 
tions fail in describing the low-energy collective modes of 
vortex configurations. These possibilities are discussed 



in more detail in Sec. IV where we discuss the shortcom- 



ings of the above-mentioned models in describing vortex 
states, and stress the need for an analysis within a self- 
consistent formalism that takes into account the dynam- 
ics of the condensate and the noncondensate on an equal 
footing. 

The structure of this paper is the following. In Sec. || 
the general formalism of time-dependent Hartree-Fock- 
Bogoliubov-Popov equations in a rotating frame of ref- 
erence is presented. In Sec. [II this formalism is applied 
to the case of a condensate containing a precessing off- 
axis vortex line, and a lower bound for the quasiparticle 
energies is derived. Furthermore, by applying this lower 
bound to the anomalous modes localized in the vortex 
core region, we show that such modes are lifted to pos- 
itive energies as the noncondensate density in the core 
exceeds a certain value. Finally, the implications and 
proper interpretation of these results are discussed in Sec. 

fvl- 



II. GENERAL FORMALISM 

We model a dilute bosonic gas with the usual effective 
grand-canonical Hamiltonian, which takes the second- 
quantized form 



-h| / drV^(r,i)V'^(r,i)V'(r,t)V'(r,t) (1) 
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in a frame rotating with the angular velocity ft w.r.t. 
an inertial coordinate system. Above, ip{r, t) denotes the 
boson field operator, and 



(2) 



is the one-particle Hamiltonian, with Vtr (r) denoting the 
external trapping potential, fi the chemical potential, and 
L the angular momentum operator. The interaction be- 
tween the particles is modelled by the effective contact 
potential Vint{r — r') = 5<5(r — r'); with the coupling con- 
stant g related to the vacuum s-wave scattering length a 
through g = ATrTi^a/m. The dynamics of the system is 
determined by the Heisenberg equation of motion 

in^il>{v,t) ^nn{Y)^{Y,t) + gi,\Y,t)^{Y,t)i,{Y,t). (3) 
ot 

Following the spontaneous symmetry breaking approach, 
we use the Bogoliubov decomposition 



V'(r,t) $(r,i) +7/'(r,i) 



(4) 



to split the field operator to the sum of the c-number 
condensate wave function <i>(r, i) — {ip{r,t)) (with (• • •) 
we signify the time-dependent nonequilibrium average) 
and the noncondensate, i.e., excitation field operator 
ip{r,t). The expectation value of Eq. (^) yields, when 
terms containing the cubic noncondensate operator prod- 
uct or the anomalous average (ipip) are neglected ac- 
cording to the Popov mean-field scheme, the generalized 
Gross-Pitaevskii (GP) equation 

d 

in— $(r,t) =£s^(r,i)$(r,i) - .9ne(r,t)$(r,t) (5) 

for the condensate wave function. Above, Cfi{r,t) = 
Ti-uijc) + 2gn{r,t), and the condensate, noncondensate 
and total gas densities are, respectively, denoted as 



ne(r,i) = |$(r,i)p, 
n{r,t) = (^t(i.,i)^(r,i)), 
n{r,t) = nc(r,t) -I- n(r,t); 



(6a) 
(6b) 
(6c) 



the chemical potential fi is implicitly determined by the 
condition 



N 



dr n{r, t). 



(7) 



where N is the total number of particles in the system. 
Furthermore, Eqs. (§)-(^ imply in the Popov scheme the 
equation of motion 

ih^^i'iv^t) = Cn{Y,t)i,{Y,t)+g<^^{Y,t)i,\v,t) (8) 

for the noncondensate operator. This equation can be 
diagonalized by expressing the field operator in terms of 



bosonic quasiparticle operators a„ and a\^. The Bogoli- 
ubov transformation 



V'(i",i) = ^K(r,i)a„ - <(r,t)aj, 



(9) 



converts the operator equation (^) to the time-dependent 
Hartree-Fock-Bogoliubov-Popov (HFBP) equations 

d 

iTi—Un{v,t) = Cn{r,t)un{v,t) -5$^(r,i)w„(r,t) (10a) 

-ih—vn{v,t) = £_n(r,t)w„(r,t) -.g$*^(r,t)w„(r,t) (10b) 

for the quasiparticle amplitudes u,i(r, i) and u„(r,t). In 
addition, for the Bogoliubov transformation to be canon- 
ical, we must impose the condition 



dr[ii*(r,t)M„'(r,t) -u*(r,t)w„'(r, i)] = (5„„' (11) 

for the amplitudes — this time-independent normalization 
can straightforwardly be shown to be consistent with Eqs. 
(pO|). Furthermore, Eqs. (^) and (^ imply for the non- 
condensate density the self-consistency relation 

= ^ [/««'("««"' + - 2Re{5„„/u„w„/} 

nn' 

+5nn'\VnW (12) 



where 



fnn'{t) = (a,\ari') , ffrm' (0 = 



(13) 



denote the normal and anomalous quasiparticle distri- 
bution functions; here and henceforth we often do not 
explicitly denote the arguments of functions, unless they 
are needed for clarity. 

Neglect of the anomalous average in the Popov scheme 
is a rather crude method to obtain a theory with gap- 
less spectrum in the homogeneous limit, as required by 
Goldstone's theorem. Effects of the anomalous average 
can be taken into account by renormalizing the two-body 
interaction potential to a T-matrix including many-body 
effects. Effectively, this can be done by replacing the 
coupling constant g by suitable spatially varying cou- 
pling functions in the mean-field equations [^41|. In 
the following analytic treatment we assume, for simplic- 
ity, the effective interaction to be constant, and only fi- 
nally briefly discuss the more general situation. 



III. PRECESSING VORTEX LINE 

Consider a repulsively interacting boson condensate 
containing a precessing vortex line, trapped by a static, 
axisymmetric potential. For simplicity, we assume the 
vortex to move in a circular orbit about the trap sym- 
metry axis, with a constant angular velocity $1 directed 
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along the trap axis. Such a precession motion has been 
observed in the experiments, and it is in agreement with 
vortex dynamics given by Eq. (||) when the nonconden- 
sate density is neglected [Q. In addition, we assume 
the system to be thermalized in the sense that the quasi- 
particle distribution functions and gnn' are constant 
in time. In such a case, the mean-field Hamiltonian is 
time-independent in the frame rotating with the angular 
velocity ft w.r.t. the laboratory coordinate system. Con- 
sequently, in such a frame Eqs. (||) and reduce to 
the stationary equations 

/:n(r)$(r) -57ic(r)$(r) =0 (14) 

for the condensate wave function and 

Cn{r)un{r) - g^\r)vn{r) = £;>„(r) (15a) 
C-n{r)vn{r) - <?$*'(r)u„(r) = -^^'"(r) (15b) 

for the quasiparticle eigenenergies E"^ and amplitudes 
Mn(r), Vn{r). Although stationary, these equations are 
not to be confused with the adiabatic approximation for 
systems with slow time dependence [^|. The adiabatic 
approximation fails if the quasiparticle states do not fol- 
low the moving vortex core rigidly, but the eigenequa- 
tions ( [Tsl ) describe the time dependence of the system in 
principle exactly. The essential difference between these 
formalisms is the terms • L contained in Eqs. (15), 
but absent in the equations corresponding to the adia- 
batic approximation — it is these terms that account for 
the deformation of the vortex structure due to its motion 

Although considerably less cumbersome than the gen- 
eral time-dependent HFBP formalism, the stationary 
Eqs. (^5|) have not been formulated nor solved so far 
for self-consistent solutions describing off-axis vortices. 
However, in the following we show that the essential ques- 
tion concerning the local energetic stability of such vortex 
states can be studied analytically under very general cir- 
cumstances. It is to be noted that the thermalization 
and the energetic stability of the system are determined 
by its excitation spectrum in the nonrotating, i.e., lab- 
oratory frame (LF) of reference, since we assume that 
the gas is confined with a static external potential and 
that the thermal component as a whole is nonrotating 
psf . Hence, in order to investigate the issue of local 
stability, we have to study the quasiparticle energies in 
the LF, instead of the rotating frame eigenenergies E'^. 
Strictly speaking, as the mean-field Hamiltonian in the 
LF has explicit time dependence via the mean fields, it 
does not possess stationary solutions, and, consequently, 
a well-defined spectrum — in the laboratory frame only 
the expectation values of the quasiparticle energies are 
well-defined. For the same reason, such a system can 
not attain exact thermodynamic equilibrium. However, 
if the thermalization rate of the gas substantially exceeds 
the kinetic rates related to the motion of the vortex, the 
system is thermalized in the sense that its dynamics be- 
comes adiabatic to a high degree of accuracy Ej, and 



the quasiparticle distribution functions approach time- 
independent values. 



A. Lower Bound for Quasiparticle Energies 

In order to investigate the positivity of the energies 
of the lowest quasiparticle excitations, we first derive a 
lower bound for the spectrum of a trapped BEC within 
the Popov approximation. 

Let ($,u, w) be a self-consistent sohition of Eqs. (|l^ ) 
and (Eq), corresponding to an eigenenergy E' and a 
given set of time-independent quasiparticle distribution 
functions /„„' and gnn'- Eqs. (|l|) and (||) imply that 
the expectation value (E) of the quasiparticle energy 
in the nonrotating laboratory frame is related to the 
eigenenergy E' in the rotating frame by the relation 
(E) = E' + il,-{V), where (1) denotes the expectation value 
of the angular momentum of the quasiparticle state. By 
expanding the expectation value / dr{tp''Ltl>) of the angu- 
lar momentum related to the noncondensate in terms of 
the quasiparticle amplitudes and distribution functions, 
one finds that the angular momentum of the system is 
changed by the amount / dr {u*Iju + vLv*) as the cor- 
responding quasiparticle occupation number is increased 
by one. Thus, we find 



E' + J dr [u*{n ■ L)u + v{n ■ L)v*]. 



(16) 



On the other hand, integration of the sum of Eq. (15a) 
multiplied by u*(r) and the complex conjugate of Eq, 
(151) multiplied by v{r) yields 



J dr [u*{£o - rj • L)u + v{Co - ft ■ L)v* - 2g<i'^u*v] 



E' / dr{\u\^~\v\^)^E', (17) 



where we have used the relation Cn = Co — ^ ■ li, the 
fact that L* = (r X fiV /i)* = — r x hV/i = -L, and the 
normalization condition (O) for the quasiparticle ampli- 
tudes. Now, combination of Eqs. (^6|) and (|l7| ) yields the 
relation 



dr {u* Cqu -\- vCqv* ~ 2g$ u*v) 



(18) 



for the quasiparticle energy in the laboratory frame. It 
is to be noted that this expression has only an implicit 
dependence on the vortex precession velocity ft, via the 
fact that the amplitudes u(r) and v{r) satisfy Eqs. (|l^). 

In order to derive a simple lower-bound expression for 
the quasiparticle energy, we first note that f dr ^'^u*v is 
real valued — this can be inferred from Eq. (|l8|), where all 
the other terms are real valued MM . Thus we have 



J dr<^^u*v < J dr \<^\^\u\\v\ 



< - I drndlu 



' + \v\') 



(19) 
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since — |u|)^ = \u\'^ + |up — 2|wHw| > 0. Combina- 
tion of this upper bound with Eq. ) yields, due to the 
positivity of the coupUng constant g, the inequaUty 

{E) > J dr [u*jCou + vCqv* - + \v\'^)g72c\ 

= j dr [u*{-h^V^/2m)u + v{-h^V^/2m)v*] 

+ [ dr{\u\^ + \v\^){Vtr+gn, + 2gfi- fi). (20) 



The kinetic energy terms may be discarded from this 
lower bound due to their positivity, and we find 



{E) > J drp {Vtr + guc + 2gfi - ^), 



(21) 



where p(r) — |u(r)|^ + |w(r)|^ denotes the density dis- 
tribution of the quasiparticle state. Using the Thomas- 
Fermi (TF) approximation 



ii^^) = ^(/i-X4r-2gn) 



(22) 



for the condensate density in the strong- interaction limit, 
we can write the second factor in the integral of the lower 
bound (|ll) as 



Vtr + gric + 2gn - fi = g{nc 



fTF) 



(23) 



Since we are mainly interested in the anomalous quasi- 
particle states which are localized in the vicinity of the 
vortex core, it is useful to split the spatial integral in 
the inequality ( pT[ ) to contributions from the core re- 
gion C and its complement C. To be specific, we de- 
fine the vortex center line to consist of points in which 
4'(r) = 0, and the core region C to contain points whose 
distance from the center line is less than the healing 
length ^ = (87rnoa)~^/^, where no denotes the local total 
density of a corresponding irrotational system. Espe- 
cially, in the strong interaction Thomas-Fermi limit we 
have 



nQ~ - {p- Vtr), 



9 



(24) 



provided that the temperature is not so close to Tbec 
that the thermal gas component would dominate over 
the condensate. Utilizing the positivity of p, Uc and g, 
we find Eqs. (21) and (p3) to imply 



{E) > / drp{Vt, + 2gfi- fi) 
Jc 

— £(max{n^'^^^ — Tic} / drp, 
C" Jc 



(25) 



where max{- • -^c' denotes the maximum value in the re- 
gion C . Furthermore, by defining the quantity 



Ic" drp 



(26) 



describing the degree of localization of the quasiparticle 
state (u, v) into the vortex core region, we can write the 
lower bound (|2^) for the quasiparticle energy in the more 
compact form 



{E) > / drp Vtr + 2gii- p, 



c 



nc}/w 



(27) 



It is to be noted that this lower bound holds for all quasi- 
particle states, when the corresponding density distribu- 
tions p and localization degrees w are used in the inte- 
grand. 



B. Anomalous Modes 

Suppose that the condensate containing a vortex line 
has an anomalous quasiparticle mode (AM) with a neg- 
ative energy in the laboratory frame. In such a case 
the system is energetically unstable, because it can, in 
the presence of dissipative mechanisms, lower its free en- 
ergy by exciting the AM with matter originating from 
the condensate. It has been argued that such "drain" 
of matter into the AM would lead to the formation of a 
binary condensate from the vortex core, and a simulta- 
neous spiralling of the vortex out of the condensate 
However, it is to be noted that when the noncondensate 
mean field is taken into account in a self-consistent man- 
ner, the quasiparticle energies depend on the distribution 
functions /„„' and gnn' ■ Especially, excitation of an AM 
increases the corresponding occupation number, and thus 
changes the energy of the mode itself. If the change in 
the energy is large enough, the anomalous mode may at- 
tain positive energies, in which case the "collapse" of the 
condensate would cease. Consequently, the vortex would 
be stabilized by the self-interaction of the noncondensate 
in the core. 

More specifically, since the quasiparticle density distri- 
bution p is positive, Eq. (^) implies that a given mode 
can be a negative-energy AM only if the expression in the 
brackets is negative somewhere in the core region, i.e.. 



< ^ (/^ - Vtr + .gnmx{n(TF) 



(28) 



somewhere in C . 

In order to estimate the vortex core localization de- 
grees w of the anomalous modes, we note that the quasi- 
particle motion is essentially determined by the potential 



off 



Vtr + 2g(nc + n) - p 



(29) 



— this is the effective Hartree potential appearing in Eqs. 
( pO| ) and (p^. For noncondensate density distributions 
not too far from the form corresponding to a thermal- 
ized state, V^s is negative in the vortex core region, and 
positive elsewhere P,|2| . Consequently, for the energy of 
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the mode to be negative, its density has to be concen- 
trated in the core region, thus yielding a locaUzation de- 
gree It; ^ 1. This heuristic reasoning for the locahzation 
degree is supported by exphcit numerical computations 
for axisymmetric vortex states of axially harmonically 
trapped condensates, resulting in values w ~ 0.7 in the 
large- limit, and w w 1.0 for systems in the noninteract- 
ing limit 1 46 1 . It is to be noted that for vortex precession 
velocities that substantially exceed the adiabaticity limit 
for rigid vortex core motion |p7| , the quasiparticles may 
be deformed and the localization degree can be smaller 
than the estimates given above. However, for the anoma- 
lous modes it always has to be of order unity. 

For harmonically trapped gases, the effective strength 
of the interactions relative to the kinetic energy in the 
system is characterized by the dimensionless quantity 
Na/ah_o, where aho = {h/miDy^^ is the harmonic os- 
cillator length, and u) = {uixi^y^^zY^^ denotes the ge- 
ometric mean of the harmonic oscillator frequencies. 
The Thomas-Fermi approximation for the condensate 
density is accurate in the interaction-dominated regime 
iVa/oho ^ 1- However, as the corrections to the chemi- 
cal potential given by the TF value ^tf = ff'^peak, where 
'^peak denotes the maximum density of the gas, are of the 
small fractional order (ISA/'a/oho)"'*^^ ln(15A^a/aho)^^^ 
its accuracy is in general on the order of 10% 
when Na/aiio ^ 1- Thus, in this parameter regime we 
may use the estimate ( p^ ) in the condition ( |28| ) . Further- 
more, as the localization degree w of the anomalous mode 
is of order unity, the contribution given by the last term 
in the inequality (|2^) is small, and we get an approximate 
condition 



(30) 



for some r G C, when Na/ai^o ^ 1 and the temperature 
is not too close to Tbec. In other words, since rio(r) is 
approximately the local maximum density in the vicinity 
of the vortex, the anomalous modes are lifted to positive 
energies as the noncondensate in the core fills half of the 
core volume. 

This result, which applies also to unharmonically 
trapped condensates, shows that within the PA even par- 
tial filling of the vortex core with non-condensate is suffi- 
cient to "suppress" the anomalous modes by lifting their 
energies to positive values. In fact the approximate in- 
equality ( |30| ) is rather pessimistic in the sense that, ac- 
cording to exact numerical solutions, even much lower 
noncondensate core densities suffice to lift the anomalous 
modes to positive energies The analysis can also 

in principle be straightforwardly generalized to cover the 
Gl and G2 approximations, if and when the effective spa- 
tially dependent coupling g is not substantially reduced 
in the vortex core region. Returning to the question of 
vortex stability, these results show that the "collapse" of 
the condensate to the anomalous core modes ceases as the 
noncondensate density due to these modes partially fills 
the vortex core and the energies of the instability modes 



rise above the energy of the condensate state. In the pres- 
ence of strong enough dissipation, this self-stabilization 
is a direct consequence of thermalization. All in all, we 
find that Popov-like self-consistent approximations pre- 
dict vortices to become locally energetically stable in the 
course of thermalization. 



IV. DISCUSSION AND CONCLUSIONS 

The Bogoliubov approximation consisting of the mean- 
field Gross-Pitaevskii equation and the Bogoliubov 
normal-mode equations has been widely used to describe 
the properties of vortices and, especially, their stability 
in dilute atomic BECs. Predictions derived from it for 
the vortex precession direction and frequency, together 
with the critical trap rotation frequencies for vortex nu- 
cleation, agree well with experimental results. However, 
the BA neglects the effects of the noncondensate. At 
temperatures T > Tbec/S a considerable fraction of the 
gas is noncondensed and the BA is insufficient to model 
the system. Moreover, we argue that for vortex states 
its accuracy is questionable even at lower temperatures. 
This is due to the tendency of the condensate matter 
to be transferred to noncondensate occupations of the 
anomalous negative-energy core modes. Consequently, 
as the system is thermalized, the vortex core becomes 
occupied by a substantial noncondensate density even in 
the zero-temperature limit. 

The Popov- like approximations, on the other hand, 
are finite-temperature mean-field theories that self- 
consistently take into account the thermal gas compo- 
nent. Their predictions for the excitation frequencies of 
irrotational condensates are in excellent agreement with 
experimental results at temperatures T Tbeo/2, but 
at higher temperatures discrepancies become obvious, as 
the noncondensate density becomes more dominant. This 
failure of the PLAs near Tbec has been attributed to the 
fact that they do not properly treat the dynamics of the 
noncondensed gas. 

On the basis of explicit numerical solutions for axisym- 
metric vortex states [OpGl and the analysis presented in 
this paper (see also Ref. |[29||, p. R177), we conclude that 
PLAs imply vortices in sufhciently properly thermalized 
and repulsively interacting condensates to be locally en- 
ergetically stable in the sense that the excitation spec- 
trum of the system is positive definite. Especially, these 
approximations imply that the vortex precession mode 
has a positive energy, which in turn seems to imply that 
off-axis vortices precess in the direction opposite to the 
vortex flow. However, this prediction contradicts exper- 
imental results, which rather show that the majority of 
vortices precesses in the direction of the condensate flow. 
In the light of the success of the mean-fleld approxima- 
tions in describing the physics of dilute condensates, the 
failure of the mean-field theory itself is an unlikely rea- 
son underlying this discrepancy. Thus we are essentially 



6 



left with two alternatives: Either the vortex states ob- 
served in the experiments have not been properly ther- 
malized, or the PLAs drastically fail in describing the 
lowest collective modes of the vortex states. In fact the 
observations of vortex precession have so far been con- 
ducted under conditions in which dissipation in the gas is 
very weak, and the degree of thermalization probably low 
1^,0. It is to be noted that, due to improper thermal- 
ization, the anomalous mode occupations may be negligi- 
ble and the vortex core almost void of noncondensate — in 
such a case the Bogoliubov approximation should natu- 
rally yield more accurate results than the PLAs. 

In any case, there is also a reason to expect PLAs to 
fail in predicting the lowest quasiparticle modes for vor- 
tex states even in the zero-temperature limit. As the 
vortex state is thermalized, the noncondensate concen- 
tration in the vortex core region becomes nonnegligible 
and, as we have seen, it substantially modifies the en- 
ergies of the lowest quasiparticle states which are essen- 
tially localized in the core. However, when describing 
the collective modes of the system, the PLAs do not 
correctly take into account the dynamics of this non- 
condensate, but treat the thermal gas peak in the core 
as a static pinning potential for the vortex. The stabi- 
lizing effect of such an "external" potential is naturally 
stronger than that of a realistic, dynamical nonconden- 
sate peak, and the lifting effect of the anomalous modes is 
consequently exaggerated within the PLAs. For the same 
reason, the prediction for the precession direction given 
by the PLAs is compromised, because the lowest mode 
rather describes precession of the vortex about a static 
noncondensate peak than the combined motion of the 
condensate and the noncondensate (see also Ref. [p9[). 
Such a failure of PLAs may be somewhat surprising, as 
one in general would expect these approximations to be 
unreliable only near Tbec, where the noncondensate frac- 
tion is substantial together with the resonant contribu- 
tions to the self-energies When vortex states are 
considered, however, the dynamics of the thermal gas 
should be taken into account at all temperatures, due to 
the filling of the vortex core with noncondensate in the 
course of thermalization. 

Consequently, the effects of the thermal gas concen- 
trated in the vortex core are not properly modelled within 
the BA nor the PLAs. In order to reliably find out the 
structure, spectrum and stability of quantized vortices 
in dilute BECs one has to use a model which treats the 
dynamics of the condensate and the thermal gas on an 
equal footing. Such fully dynamic approximations for in- 
homogeneous gases have been developed in recent years 
p9|-p4| , but they have not yet been applied to study vor- 
tex states. Also due to the advances in experimental tech- 
niques and measurement accuracies there is a need for 
investigating vortex states beyond the Bogoliubov level. 
Because of the strong inhomogeneities, vortices are po- 
tentially quite efficient testbenchcs for thermal field the- 
ories in modelling atomic BECs. 

In conclusion, we have studied the stability of precess- 



ing off-axis vortices in dilute, trapped Bose-Einstein con- 
densates within the self-consistent Popov approximation 
and its extensions Gl and G2. By deriving a lower bound 
for the excitation energies, it has been shown that within 
these approximations the vortex states are locally ener- 
getically stabilized as the vortex core is partially filled 
with thermal gas. Furthermore, as such filling of the vor- 
tex core with noncondensate is an inherent consequence 
of thermalization and the existence of anomalous modes 
for vortex states in pure condensates, we argue that the 
effects of the thermal gas for the physics of vortices may 
not be neglected even in the zero-temperature limit. Fur- 
thermore, the shortcomings of the Bogoliubov approxi- 
mation and the Popov-like self-consistent approximations 
in modelling vortex states have been pointed out and 
analysed. Essentially, the BA may be too "pessimistic" 
in predicting vortices to be definitely unstable, and the 
PLAs too "optimistic" in predicting them to be self- 
stabilized by the thermal gas component. We argue that 
a reliable analysis of the stability of quantized vortices 
in dilute BECs requires a fully dynamic formalism — such 
an investigation remains a challenge for further studies. 
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